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Abstract 

Consider a three dimensional system which looks like a cross-connected pipe system, i.e. a 
small sample coupled to a finite number of leads. We investigate the current running through 
this system, in the linear response regime, when we adiabatically turn on an electrical bias 
between leads. The main technical tool is the use of a finite volume regularization, which 
allows us to define the current coming out of a lead as the time derivative of its charge. We 
finally prove that in virtually all physically interesting situations, the conductivity tensor is 
given by a Landauer-Biittiker type formula. 



1 Introduction 

In any experiment which involves running a current through a microscopic sample, having a good 
quantum description of various conductivity coefficients is a priority. Such a description has been 
first derived by Landauer [5T] and then generalized by Biittiker [3]; this is what one now calls 
the Landauer-Biittiker formalism. The main idea is that the conductivity of mesoscopic samples 
connected to ideal leads where the carriers are quasi free fcrmions, is completely characterized by 
a one particle scattering matrix. 

Many people have since contributed to the justification of this formalism, starting from the 
first principles of non-equilibrium quantum statistical mechanics. In this respect, there are two 
different ways to model such a conduction problem, and let us briefly discuss both of them. 

The first approach is the one in which one starts with several decoupled semi-infinite leads, each 
of them being at equilibrium [19j . Let us assume for simplicity that they are in grand canonical 
Gibbs states having the same temperature but different chemical potentials. Then at t = they 
are suddenly joint together with a sample, and the new composed system is allowed to evolve freely 
until it reaches a steady state at t = oo. Then one can define a current as the (Cesaro limit of the) 
time derivative of a regularized charge operator, and after lifting this regularization one obtains 
the L-B formula. This procedure is by now very well understood and completely solved in a series 
of very recent works (see pQ, |25j and references therein). Let us mention that in this approach 
a crucial ingredient is the fact that the perturbation introduced by the sample and coupling is 
localized in space. One can even allow the carriers to interact in the sample |18| . and the theory 
still works (even though the L-B formula must be replaced by something more complicated). Note 
that even if we choose to adiabatically turn on the coupling between the semi-infinite leads, the 
result will be the same. A rigorous proof of this fact is in preparation |14j . 
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The second approach is the one in which the leads (long, but finite) are already coupled with the 
sample, and at t = — oo the full system is in a Gibbs equilibrium state at a given temperature and 
chemical potential. Then wc adiabatically turn on a potential bias between the leads, modeling in 
this way a gradual appearance of a difference in the chemical potentials. The statistical density 
matrix is found as the solution of a quantum Liouville equation. The current coming out of a 
given lead is defined to be the time derivative at t = of its mean charge. Then one performs 
the linear response approximation with respect to the bias thus obtaining a Kubo-like formula 
[6], and finally the thermodynamic and adiabatic limits. The current is given by the same L-B 
formula, specialized to the linear response case. Note that in contrast to the previously discussed 
approach, the perturbation introduced by the electrical bias is not spatially localized, and this 
makes the adiabatic limit for the full state (i.e. without the linear response approximation) rather 
difficult. 

Some significant papers from the physics literature which initiated the second approach are 
[TU] . [2"T] . [2"8] , [T7], [H] and [5]- They contain many nice and interesting physical ideas, even for 
systems which allow local self-interactions, but with no real mathematical substance, also due to 
the fact that many techniques were not yet available at this time. 

A first mathematically sound derivation of the L-B formula on a discrete model was obtained 
in [12j and further investigated in [13] . In the current paper we greatly improve the method of 
proof of [Hj, which also allows us to extend the results to the continuous case. 

There are many interesting and hard questions which remain to be answered. Namely, can one 
compute the adiabatic limit in the second approach without the linear response approximation? 
If yes, then what is the connection with the first approach? Are the two steady states identical? If 
not, on which class of observables do they have equal expectations? Can one say anything rigorous 
about transient currents (see [23] and references therein)? 

A distinct issue is the time dependent coupling and/or bias introduction, eventually periodic in 
time, closely related to the pumping problem. Wc here mainly refer to the so called BPT formula 
[H], rigorously investigated and justified in a series of papers by Avron et al [2J, [3]. We also note 
that Landaucr type formulas arc universal in a way, and can be met in various other contexts like 
for example in the viscosity experienced by a topological swimmer [4]. 

The paper is organized as follows. 

Section 2 introduces the mathematical model, and derives a current formula (|2.44[) via linear 
response theory, finite volume regularization and adiabatic switching of the electrical bias between 
leads. 

Section 3 contains the proof of the thermodynamic limit (i.e. the length of the leads is taken 
to infinity). The new current formula is given in (|3.25p . 

Section 4 contains the adiabatic limit of the current, and makes the connection with the 
stationary scattering objects in (|4.14[) . 

In Section 5 we derive the usual Landauer-Biittiker formula from (|4.14p : we state our main 
results in Theorem 15.21 Wc also derive a continuity equation in subsection 5.1, which shows that 
the steady current is the same, no matter where we measure it on a given lead. 

2 The finite volume regularization 

Since our results can be easily extended to one and two dimensions, wc will only consider the three 
dimensional case. The small device through which the current flows is modelled by a bounded 
domain of R 3 , not necessarily simply connected, but with a regular enough boundary. This is the 
" sample" , which is linked to a finite number of " leads" (which can be cylinders of length L and 
radius 1). Let us for simplicity only consider two leads, both cylinders being parallel to the first 
coordinate axis and located in — L — a < x\ < —a (the left lead) and a < x\ < L + a (the right 
lead). Here a is a nonnegative number. The sample is a subset of the slab —a < x\ < a. Note 
that one possible system is just the union of the two leads, when a = and the sample is absent. 
If a > we demand the sample to be smoothly pasted to the leads, so that the boundary of their 
union is at least C 2 . 
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2.1 The equilibrium state 

The union between the leads and sample will be denoted by Xl- The one particle Hamilton 
operator is given by Hi := — Al +w, where — Al is the Laplace operator with Dirichlet boundary 
conditions in Xl- The potential w is smooth and compactly supported in the region where the 
sample is located, i.e. Xl fl {x £ I 3 : —a < x\ < a}. Without loss of generality, we put w > 0. 
The one particle Hilbert space is TL L := L 2 (Xl). 

Denote by M D (X L ) := H£(X L ) n H 2 (X L ), where H&(X L ) and H 2 (X L ) are the usual Sobolev 
spaces on the open domain Xl C R 3 . Since we assumed enough regularity on the boundary of 
Xl, the operators: 

- A L :M d (X l ) H l (2.1) 
H L := -A L + w(Q) : M D (X L ) -» H L (2.2) 

are self-adjoint on Md(Xl). It is also very well known that their spectrum is purely discrete and 
accumulates to +oo. Due to our assumptions on w, the Hamiltonian in (|2.2|) is a positive 
operator and its resolvent is denoted by 

R L (z) := (H L - z)~ Y (2.3) 

for all z € C \ R+; we shall simply denote R^ := Rl(— 1). 

We only consider the grand-canonical ensemble. In the remote past, t — > — oo, the electron gas 
is in thermodynamic equilibrium at a temperature T = \ > and a chemical potential fi. The 
appropriate framework is the second quantization. Here the Hilbert space is $l '■= © (Tt L ) An , 

nGN 

where (H L ) An is the n times antisymmetric tensor product of TL L with itself. Let = dT(li) 
be the number operator (where 1l is the identity operator on TL L ). 
The second quantization of H l is denoted by \j L := dT(Hi). Let 

u L (t) := T(e-* tHL ) = e~ lti >i- 

be the associated time-evolution of the system. Under our conditions, the operator e~" UL € 
V>i[H L ], the ideal of trace-class operators on H L , and this property extends itself to the second 
quantized operators. 

The associated Gibbs equilibrium state is: 

e -0(<)i,-f«n£) 

Pl ;= Tr{e-^-^)}' (2 ' 4) 

It is a standard fact (see Proposition 5.2.23 in [7]) that if we denote by p(X) := [e 13 ^^ + l) , 
then for any bounded operator T £ B[H L ] we have: 

Tr 5i {p L dr(T)} - Tr w , {p(H L ) T} . (2.5) 

2.2 Turning on the bias 

Let us introduce three projections defined by natural restrictions on the left lead, sample, and 
right lead respectively: 

n_ : H L -> L 2 (X L n {x e K 3 : -L - a < x x < -a}), 
n+ : H L -> £ 2 (<Y L n {x e R 3 : a < .T! < L + a}), 

n : H L -> L 2 (X L n{x£l 3 : -a < xi < a}). (2.6) 
If u± £ R, define the bias between leads as 

V := v-n_ +v+n+. (2.7) 
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We will now consider the Hamiltonian describing the evolution under an adiabatic introduction 
of the electric bias. Consider a smooth switch-on function \ which fulfills XiX' G °°?0)), 
x(0) = 1 (note that the first two conditions imply linit^-oo \{t) = 0). Let r] > be the adiabatic 
parameter, and define x ?? (i) := xivfyi V n (t,x) :— Xr t (t)V(x). At the one-body level we have 

K LtV (t) := H L + V v (t, Q) = H L + x(vt)V(Q) (2.8) 

where V(Q) denotes the bounded self-adjoint operator of multiplication with the step function 
V(x). We shall use the notation Kj,^ := Kj, ! ^(0) = + V(Q). The evolution defined by this 
time-dependent Hamiltonian is described by a unitary propagator W(t, s), strong solution on the 
domain of of the following Cauchy problem: 

-id t W L ,ri(t, s) = -K L , v (t)WL,ri(t, s) , > 

W L , v (s,s) = l [Z ^> 
for (s,t) E M 2 . 

Remark 2.1. For any L < oo and rj > 0, the operators {Ki^(t)} t£ R are self-adjoint in H L , 
having a common domain equal to Md(Xl). They are strongly diffcrcntiable with respect to t € M. 
with a bounded self-adjoint derivative 

dtK L „(t) = vx(vt)V(Q). 

The following result is standard and we state it without proof: 

Proposition 2.2. The problem (|2.9|) has a unique solution which is unitary and leaves the domain 
Md(Xl) invariant for any pair (s,t) £ R 2 . For any triple (s,r,t) 6 R 3 it satisfies the relation 
WL jTI (t, f)Vl / i,.^(r, s) = WL,ri{t, s). Moreover it also satisfies the equation 

- id s W L , n (t, s) = W L , v (t, s)K LtV (s). (2.10) 

Later on we will need the following simple but important result: 

Proposition 2.3. Let us define the commutator 

[H L , W L ,r,(t, s)) = U L W L . v (t, s) - W L . v (t, s)H L 

on the dense domain Wr)(Xi,) (left invariant by the evolution). Then the commutator can be 
extended by continuity to a bounded operator on TL L such that: 

sup ||[H L) Wl,„(M)]||<oo, (2.11) 

( s ,t)eK= 

sup ||(H L + l)H/ L ^(<,s)(H L + l)- 1 || <oo. (2.12) 

(s,t)£Ri 

Proof. Let us fix some <j) m ^d{Xl), and for t < let us consider the application 

(-oo, t) 3 s i ► 0(a) := \K LlTI (s),W L , v (8,t)]<i> £ L 2 (X L ), 

K Llf ,(s)<t>(s) CU D (X L )*, (2.13) 

with (f>(t) = 0. Due to these properties, the function <fi(s) is diffcrcntiable in the weak topology on 
IId(Ax)* and we can compute its derivative: 

d a </>(s) = [d s K L ,n(s),WL, n (s,t)}(t>+ [K L ,n(s),d s W L ,n(s,t)}4> 

= TO'M[V(Q) ) ^(s,t)]0-iK £>7) ( S )[K Lj7) ( S ),W Lj7) (s ) t)]0 

= -iK Ltn (s)(f>(s) + i>(s), (2.14) 



4 



where ip(s) := 'HX f (' l l s )[^ r (Q)i Wl,t){s, t)\4>. Thus, taking into account the above solution of the 
Cauchy problem (12. 9| . we can write the solution for <p(s) as 

cj)(s) = J W L ,r,(s, r)ijj{r) dr. 

Let us recall that WL, v (t, r) is unitary, the function \' belongs to L 1 (M_) and || V"(Q)|| < max{f_, «+} 
so that 

\\<Ks)\\ h l < ||x'|| L i (K _)max{i;_,«+}||0||, Vs G (-00, t), Vt G »_. 
Thus for any pair of functions ip and </> in Md(Xl), dense in Ti. L , we have 

\(il>,\K L , v (a),WL, v (B,t)]<t>)\ < ||x'|U MR _)max{«_,« + }||V||||^||. 

After a straightforward density argument, the use of Riesz' representations theorem, and the fact 
that V(Q) is bounded, we obtain a uniform estimate as in (|2.1ip if —00 < s < t < 0. To obtain 
the same for — 00 < t < s < we remind that Wl, v (s, t) = WL. n (t, s)* and for -0 and <j> in Wd(Xl) 
we can write (using the invariance of the domain Wd(Xl) under Wi^t, s)) 

-([Hi.^tMllf^) = (0, [U L ,W L ,r,{t,s)]i>), 

and we easily get the desired conclusion. Finally, (|2.12[) is an easy consequence of (|2.11l) . □ 



Now let us consider the second quantization of the above time-dependent evolution. 

dT{K L ,, 1 {t)) = i )L + x v {t)dY{V{Q)). 

Let us first observe that the perturbation dT(V(Q)) is no longer bounded on It is rather easy 
to verify that the family {tOL iTI (t, s)}( s ,t)eR2, that is well defined by ttJj, j77 (t, s) ■= T(WL,, 1 {t, s)) 
gives the unique (and unitary) solution of the Cauchy problem: 

-id t to L>ri (t, s) = -dT (K L:V (t)) to LiV (t, s) . . 

V0L, n {s,s) = 1. [ ' ' 



2.3 The nonequilibrium state 

In our adiabatic approach, the nonequilibrium state is completely characterized by a density matrix 
which solves the Liouville equation. Its initial condition at t = — 00 is the Gibbs equilibrium state 
pL, associated to the Hamiltonian fj L at temperature and chemical potential /i. We let the 
state evolve with the adiabatic time-dependent Hamiltonian Ki jT) (t) up to t = 0. Then we perform 
the thermodynamic limit L — > 00, and finally we should take the adiabatic limit r\ — ■> 0. This 
last step raises for the moment some important technical difficulties. We will avoid them by only 
considering the linear response behavior, i.e. the linear contribution in the bias v+ — V—. The 
proof of the existence of the thermodynamic and adiabatic limit for the full density matrix is a 
challenging open problem. 

Now let us formulate and solve the Liouville equation. We first want to find a weak solution 
to the equation: 

f id t p L ,r,(t) = IdT (K L!V (t)) , p £)I) (t)], t < 0, 
\ limt^-oo p L>v (t) =p L . 

Let us note that for a fixed s, the operator (see also (|2.15|) ) 

tt>L,ri(t,s)pLtt> L>rl (t,s)* 

solves the differential equation, but does not obey the initial condition. What we do next is to 
take care of this. 
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The state pi being defined in (|2.4p commutes with any function of t) L . Thus we can write: 
m LtV (t, s)p L tt> LiTI {t, s)* = m L , v (t, s)e^- s ^^p L e-^- s ^^w L!rl {t, a)* 
= r (w L , n {t,s)e^- s ^) p L T (w LtV (t, s)e-^- s ^y, (2.17) 

and we can consider the two-parameter family of unitaries: 

Sl L , v (t,s) := W L ^t,s)e^- s)li \ (M) e R 2 . (2.18) 

We will now show that the solution at time t is given by the following strong limit (norm limit on 
each sector of fixed number of particles): 

p L , v (t) = lim T(Q Ltr ,(t,s))p L T(n L , n (t,s))*. (2.19) 

s — > — oo 

Let us first show that the limit exists. Due to the continuity of the application T, it is enough to 
prove the existence of the norm limit lim CIl n(t, s) in the one particle sector. For that we use 

the following differential equality, valid in strong sense on the domain of Hi: 

-id a (l L , v (jt, s) = W L , v (t, s)(K L>v (s) - n L )e^ s ^ L . (2.20) 
If we denote V(r) := e lrHi V(Q)e~ lrIlL we observe that 

/ -id s Q L , v (t, s) = x(v s )^L, v (t, s)V(s - t) , , 

and thus 

rt 



n L . n (t, s) = 1 - i J x(vr)n L , v (t, r)V(r - t)dr. (2.22) 

Since f2i >)7 (i, r) is unitary for any pair (s, t) G K_, \\V(r — t)\\ is uniformly bounded in r, and \ 1S 
intcgrable, we conclude that the limit 



lim Q, L , n (t,s) =: Sl LtV (t) (2.23) 

s— > — oo 

exists in the norm topology of B[7i L ] and thus defines a unitary operator. 

Proposition 2.4. The operator f2i i?? (i) preserves the domain o/Hi, and the densely defined 
commutator [Hi, f2i i?J (i)] can be extended to a bounded operator on 7i L , with a norm which is 
uniformly bounded in t < 0. 

Proof. If ip belongs to the domain of Hi, then for every <fi £ Tt L we have: 

KH L ^n L)tJ (t)(H £ + i)-V>l 

= lim |(H L ^^(t, S )(Hi + l)-V(*- s ) H ^)|<C7||V|| \\<t>\\, (2.24) 



s — > — oo 



where we used (|2.18|) . (|2.23|) and (|2.12p . Now since Hi is self-adjoint, it means that f2i i?) (i)(Hi + 
l) _1 (/> belongs to the domain of Hi, hence: 

sup ||Hitti.„(t)(Hi + l)" 1 !! < C. (2.25) 

t<o 

Now for ij) and <j) in the domain of Hi we can write 

[Hi,Oi,,(f)]0)| 

- lim |^,[Hi,M/i, ) (i, S )]e 4 (*- s)H ^)| < C|H| Ml, (2.26) 

s — > — oo 

where we used (|2.11j) . A density argument finishes the proof. □ 
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Up to now we have shown that pL,r t is a weak solution to the Liouville equation (|2.16|) . This 
density matrix is a trace class operator. The key identity which gives the expectation of a one-body 
bounded observable T G M(TC L ) lifted to the Fock space is (see also (|2.5[) ): 

Tr dL ( H , v (t)dT(T)) = Tr Si (p L T(n L , n (t))*dr(T)T(n Ltn (t))) 

= Try, { PL dT{ni v (t)Tn L , ri (t)}) 

= Tr nL (Q L , v (t)p(U L )Ql >ri (t)T) . (2.27) 

We stress that all this works because Xl is finite. The main conclusion is that if we are only 
interested in expectations of one body obscrvablcs, the effective one-particle density matrix is: 

PL At) ■= ^L,r,{t) P {U L )9,* L ^t) G Bi(H L ). (2.28) 

We can now prove that the above mapping is diffcrentiablc with respect to t in the trace norm 
topology. We write: 



PL, v (t)= lim ^L,r,(t,s)p{H L )ni(t, S ) = lim W L , rl {t, S )p(U L )W*(t,s) 

s—> — oc ' 1 s—> — oc ' ' 

= W L ,r,(t,0) ( lim W L , 1 (0 7 s)e-^p(n L )e^Wl n (0 7 s)\w L , 1 (t 7 0r 

= w L , v (t, o)pL, v (°)Wi, v (t, oy = w L>v (t, o)n LiV (o)p(n L )niAo)wi v (t, o). (2.29) 

We have Hi ( o(Hi)Hi G Bi[7Y L ] because of the exponential decay of p. Using (|2.25p we obtain 
that Hl/9l iJ7 (0)H£ is also trace class. Then using ()2.9[) and ()2.12p we conclude that 1 1— > pL,n(t) S 



li \H L ] is diffcrentiablc and: 



dtPL, n {t)\ t=0 = -i \K L , v ,p L ,r,] G Bi[H L ]. (2.30) 

As we have said at the beginning of this subsection, studying the limits L — > oo and r/ — > 
on the above formula for the whole state seems to be rather difficult, and we will only consider 
the first order correction with respect to the potential bias, obtained by considering the equation 
([2~2"2"]l in the limit s —> -oo: 

SW*) = 1 - * / x(vr)^L, v (t, r)V(r - t) dr 



— oo 

t 



X(??r)y(r - 1) dr + 0{V 2 ). (2.31) 

Let us point out here that a control of the above rest 0(V 2 ) is a difficult task that we will not 
consider for the moment. 

Having in mind the above argument, we define as the linear response state at time 0: 

p L , v :=p(U L )-[V v ,p(U L )}, (2.32) 

where: 

V v := i f xMV(r) dr. (2.33) 



2.4 The current 

The main advantage of our approach is that we can define the current coming out of a lead as the 
time derivative of its charge. We define the charge operators at finite volume, to be the second 
quantization of projections H± (see (|2.6[1 ): 

Q ± :=dT(II ± ). (2.34) 
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The average charge at time t is given by: 

q(t) := Tr 5i (|>L,„(t)0±) = Tr w , (pi,„(i)II ± ) . (2.35) 

By differentiating with respect to t and using the conclusion of the previous subsection we obtain 
the average current at time t — 0: 

iL.r, = -iTr H L ([K L!V , p L ^\ n±) 

= -iTr nL (\H l ,Pl, v ] n±) - iTr HL ([V(Q), p Li „] n±) . (2.36) 

But (see ([577])) Tr w z. ([V"(Q), pi,,,] n±) = -Tt h l {p L „ [V(Q), U±\) = and we deduce that: 

jh.-q = -iTr H L ([H L ,p L . n ] U±) 

= -iTr nL ((H L + 1K,„(H £ + l)R L n± - Rl(H l + 1)pl,„(H l + i)n±) 

= -iTr HL ((H L + 1)pl,„(H l + 1) [R L , II ± ]) . (2.37) 

Because we only are interested in the linear response, we will use (|2.32p . Remember that p(Hi)H^ G 
Bi[7Y L ]. An important observation is that the following commutator defined as a sesquilinear form 
on Wo{Xl) 2 , can be extended to a bounded operator on L 2 (Xl) since: 

[H L ,VJ = / xM {5 s e^l/e— Hi } 

= V-r) X '(vs)e lsllL Ve- lsllL ds. (2.38) 

•/ — oo 

Note that we do not have to commute with V in order to get this result. In fact [Hi,, V] is 
quite singular due to the sharp characteristic functions from the definition of V. 

Thus the second term in (|2.32j) is also trace-class and we can write the linear response average 
current at time t = as: 

j L ,r, ■= -iTr H L ((Hi + l)p L , n (H L + 1) [H L ,Il ± }) 

= -iTv nL ((Hi + i)p(n L )(n L + 1) [R L ,n±]) 

+ iTr^((H i + l)[V, 7 , i o(Hi)](H i + l) [R L ,n±]). (2.39) 

The first term of the last sum is zero due to trace cyclicity. We then obtain: 

j L ,v = ((H £ + 1) [V n ,p(H L )} (Hi + 1) [Ri,n±]) 

= iTv nL ((Hi + 1) [V v ,n L ] (H L + 1) [p(H L ),n±]) , (2.40) 

where the second equality is obtained by carefully developing the commutators and using the 
cyclicity property of the trace and the fact that H|p(Hi) is trace class. 

Now we rewrite (|2.40[) in a more suitable form for taking the limits L — > oo and 77 — > 0. 

We begin by computing the first commutator in (|2.40p : 

[V„,Rl] = * / Xiw) \v{r),H L ] dr 

J — OO 

r° 

= i xMe"' Hl [V(Q),R L }e-' lrUL dr. (2.41) 

J — OO 

Now let us observe that in the strong topology on B[WzJ we can write: 

d s e tsUL R L (z)V(Q)R L (z)e~ ts H - = ie^ [V(Q), Ri(z)] e~' sHi , 
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so that we get: 

,o 

[V V ,R L } = / xM {a se isHi R L l/(Q)R L e- isHi } ds 

J — oo 

r° 

= R L V(Q)R L - V x '(ns)e lsIlL R L V(Q)R L e~* sUL ds. (2.42) 

J — OO 

Inserting the first term in (|2.40[) , we observe that it gives zero: 

Tr nL (V(Q)[p(H L ),IL±]) = 0, (2.43) 

due to the trace cyclicity and the fact that V(Q) = v + H + + V-IL- commutes with H±. Thus 
the final expression is: 

~3L, n = -iv[ x(r 1 s) r rT HL (e lsIi -V(Q)e-^[p(n L ),n ± })d S . (2.44) 

J — OO 

3 The thermodynamic limit 
3.1 The limit of the dynamics 

For different values of L, the Hamiltonians Hi are densely defined in different Hilbert spaces 
T-l L := L 2 (Xl). In order to study the behavior of our system at L — > oo, we embed Tt L in the unique 
Hilbert space TL := L 2 (<Y 00 ) (i.e. with infinitely long leads), using the natural decomposition: 

TL = H L ®L 2 (X ao \X L ). (3.1) 

We denote by the orthogonal projection corresponding to the restriction to X^. Its orthogonal 
corresponds to the restriction to X x \ Xl. Each operator Hi is self-adjoint on the domain 
W-d(Xl) (see (JHU)), while the resolvent Rl(z) is a bounded operator. Wc extend Rl(z) to TL 
by the natural formula Hj 1 Rt 1 {z)T\i j ] this new z-dependent operator family is a pseudoresolvent, 
which vanishes on Hj^H. 
Let us denote by 

H:=-A + w(Q) (3.2) 

with —A minus the usual Laplacian with Dirichlet boundary conditions defined on the Sobolev 
space iJ 2 (<Y 00 ). We denote by R(z) the resolvent (H — z) _1 (and with R := R(— 1)). 

Proposition 3.1. The sequence of operators {niRini}i>i converges strongly to R. 

Proof. Fix u £ H 2 (Xoo). If t[) e Cq°(K), define ^(x) := tp(xi). If L is large enough such 
that supp(V>3u) C Xl, then tp^u £ Wd(Xl). If we denote by D\ the subset of ff 2 (A' 00 ) with 
compact support in the x\ variable, then H is essentially self adjoint on it. Then the subspace 
£ := (H + 1)D\ is dense in the Hilbert space TL. For any u € £ , there exists L u > such that for 
every L > L u we have IliRii = Rw, Ru <E Wd(Xl) and HiRw = HRu, and thus we have: 

n L R L n LU Ru = n L R L (n L H - H L ) Ru = U L R L Uiu. 

As ||R|| < 1 and ||Rl|| < 1 for any L > 0, the proof follows after a density argument. □ 

Corollary 3.2. The sequence of operators {Hie** Hl -ni}i >a converges in the strong topology to 
e ItH , uniformly for t in any compact subset ofM.. 

Proof. Let us fix u G TL and some We can take u with compact support, such that IliU = « 

for L large enough. For any continuous function which vanishes at infinity ip £ Coo(R) consider 
V->(Hi) defined on H L through functional calculus. We define the operator ^(Hi) = Uli/j^^IIl. 
Let us observe that Coo(M) is the norm closure of the algebra defined by the 'resolvent functions' 
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{v z (t) := (i — 2) }zgc\]R and our definition of iP(H.l) agrees with the definition of the extension 
of the resolvent R^,. Thus the strong convergence of the pseudoresolvents for L — > oo immediately 
implies the strong convergence of the operators %j}{R£) to f/>(H) defined by the usual functional 
calculus on H. Now let us choose ip G Coo(R) such that i/j(X) = 1 for |A| < 1 and ■tp(X) = for 
|A| > 2. For any E > let us denote by V'e(A) := i/j(E~ 1 X). Using the above arguments, let us fix 
e > as small as we want and choose E > such that ||(1 — -0b(H))m|| < e/2. Then there exists 
L € large enough such that ||(1 — ^(Hl))u|| < e for any L > L e . Then 

\\(Tl L e uu -Tl L ~ e UH )u\\ 

< \\(e^f E (H L ) - e^f E (H))u\\ + ||(1 - ^e(H l ))u\\ + ||(1 - <MH)H|. 

For the first term we observe that for any fixed fgR the function A i— > i/>(A)e ltA belongs to Cq(R) 
and thus we can use once again the strong resolvent convergence in order to control it. Due to 
the continuity of the map t i— > e ttx (at fixed A), we deduce that we have the strong convergence 
locally uniformly in t G R. □ 



3.2 The limit of the current 

Let us fix some small enough r > and consider the positively oriented contour: 

C r = {X + ir | A g M+} U {re 10 \ 6 G [tt/2, 3vr/2]} U {A - ir \ A G R+} C C. (3.3) 

Then by analytic functional calculus we can write (due to the analyticity and the decay properties 
of the function p): 

p(U L ) = JL J p(z)K L (z)dz. (3.4) 

All these operators can be extended to the whole Hilbcrt space by the procedure Hlp(H.l)^l, 
(i.e. when considered in 7i, p(Hl) stands for n^^H^ni) but for notational simplicity we drop 
the cut-off projectors. The main result of this subsection is: 

Proposition 3.3. The operators [p(Hi),IIJ, 1 < L < oo, are trace class. Moreover, for every 
n > 1 there exists C > such that for L > 1 we have: 

||[p(H),n+] - [ P (H L ),n + ]|| Bi(w) < CL~ n . 

Proof. Let < x( x i) < 1 be a C°° cut-off function such that x = 1 on xi > a + 2 and \ — on 
x\ < a + 1. We have n + x = Xi an d the following identity holds: 

[R L (z),u+] = [R L (z), x ] + R L (z)n + (i-x)-n + (i-x)RLM. (3.5) 

Now choose two functions $,$ G C°°(M 3 ) such that < $,$ < 1, $(x) = 1 if < 1/4, $(x) =0 
if |xi| > 1/2, and $ 2 (x) + $ 2 (x) = 1 on R 3 . Denote by $ L (x) := $(x/£) and $ L (x) := $(x/L). 
We have: 

$|(x) + $l(x) = l, VxeR 3 , (3.6) 
supp($ L ) C {|x| < L/2], supp($ L ) C {|x| > L/4}, (3.7) 

||Z? Q $ L |(x) + |£> Q $ L |(x)} < C(a)L- |Q| . (3.8) 



sup 

x£l 3 



If z G C \ R define: 

S L (z) := $ L R(z)$ L + $ L R L (z)$ L . (3.9) 
The range of this operator is included in the domain of Hz, and we can write: 

(H L -z)S L (z) = l + T L (z), (3.10) 

T L (z) = [-2(V$ L )V - (A$ L )]R(z)$ L + [-2(V$ L )V - (A$ L )]R L (z)l> L . (3.11) 
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This leads to a resolvent-like equation: 

R L (z) = S L (z)-R L (z)T L (z). (3.12) 

Using standard Combes-Thomas estimates |llj . [15] . one can prove the following lemma (given 
without proof): 



Lemma 3.4. Denote by d(x) := ^/l + x\ . If z € C r (see (|3.3|) ) we denote by X = + 'St(z) 2 . 
Fix < /?2 < 0i < 1. Then there exist 8 > 0, p > and C > such that uniformly in L > 1 
(including L = oo) ) and z 6 C r : 

\\e±™D a R L {z)e^\\ B{m + ||e ± *Ri(«)e = F*|| B(ia , i o 0) < CA*, |a| < 1, (3.13) 
||e-^"R L (z)e & "|| BS <CAf. (3.14) 

We are now ready to prove Proposition 13.31 The idea is to use ()3.5p in (|3.4p and prove that 
each term on the right hand side converge in the trace norm to the operator with L = oo. 

Let us first prove that p(H)II+(l — x) G Mi(H). Because p has an exponential decay at infinity, 
we can integrate arbitrary many times by parts in (|3.4p and replace R{z) by R w (z), N > 2: 

p(U) = JL J p N (z)R N (z)dz, (3.15) 

where p^ still decays exponentially at infinity. Take N = 2. Let us prove that the integrand in 
(|3.15p becomes trace class if we multiply it by the localized multiplication operator <fi := 11+ (1— %). 
Indeed, for < ct2 < ol\ < 1 we have: 

R 2 (z)0 = {R(z)e~ Q2 ^} { e Q ^R(z)e- Q1 ^} e Ql ^^. (3.16) 

The operator e ai ^(f> is bounded, while the other two are Hilbert-Schmidt (see Q3.14P ). Thus: 

\\R 2 (z)4>\\ Bi{h) <C\p, VzeC r . (3.17) 

Now since pn has an exponential decay, the integral with respect to z defines a trace class operator. 
Let us now prove that for any n > 1 there exists C > such that: 

||n+(i - x )p(H) - n + (i - x)p(Kl)\\ MH ) < c L ~ n - (3-is) 

Using $iIT + (l — x) = and $ill + (l — x) = n + (l — \) for L > Lq, and introducing equations 
((3~1"2"|) and (|33)) in (|3~4) we obtain (denote again by = 11+ (1 - x))- 

0p(H)-0p(H L ) = 0p(H)($ L -l) + ^ / p{z)4>R L {z)T L {z)dz. (3.19) 

Both operators on the right hand side are trace class with a fast decaying trace norm. For the 
first one we use ()3.15p with N = 2 and reason as in (|3.16p : 

0R 2 (z)(l-$ L ) (3.20) 
= e Ql ^{ e - Q ^R(z)e Q ^} {e- a ^R(z)e a ^} e^^ {1 - $ L ), 
< otz < a.2 < ai < 1. 

Again e ai ^~(f> is bounded, the next two operators are Hilbert-Schmidt, and the last one has a norm 
bounded from above by e~^, c > 0, because of the support properties of 1 — Hence 

||0R 2 W(l-^)|| Bi(R3) <CA^-^ 
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and after integration with respect to z we obtain a decay as in (|3.18|) (use the exponential decay 
of pn)- 

As for the other term in (|3.19p . we have to integrate by parts with respect to z in order 
to obtain products containing three resolvents. One resolvent is used to bound the momentum 
operator from the formula of Tl, and the other two for creating two Hilbcrt-Schmidt factors. Let 
us consider one typical term: 

0R L (z)$' L 9 Xl R|(z). 

We write: 

^R L (z)^ L d Xl Rl(z) = {e Ql ^} {e- Ql "R L (z)e Q2 "} 
. {e- a *%& L } {e- a *%d xl R L (z)e a *%}e- a ^ll L (z), 

< a 2 < ai < 1. (3.21) 

The first factor in the above product, e Ql ^</>, is bounded. The second and the fifth factors are 
Hilbcrt-Schmidt, while the fourth one is bounded (see again (|3 . 1 3|) and (|3.14p ). Now the operator 
e~ Q2 2X$' L i s again bounded by since <&' L is supported in \x\\ > L/A. Hence the operator in 

(|3.2ip has an exponentially small trace norm, and after integration with respect to z we obtain 
an estimate as in the right hand side of ()3.18j) . 

Looking back at (|3.4p and (|3.5p we see that we also need to prove that [p(H), x] is trace class, 
and 

||[p(H), X ] - [p(H L ), X ]\\ Vi(n) < CL- n . (3.22) 

Since (use (|H7T5|> with N = 2) 

[p(U), X } = ^-[ p 2 (z){R 2 (z)[x,H]R(z) + R(z)[ X ,H]R 2 (z)} dz (3.23) 

111 JC r 

we see that a typical contribution to the integrand is an operator of the form H(z)x'd xl H 2 (z). 
We can write (0 < a\ < 1): 

n(z)x'd Xl n 2 (z) (3.24) 

= {H(z)e- a ^} {x'e 2a ^} {e- a ^d Xl R(z)e ai ^} {e^^R(z)} 

The first and last factors arc Hilbcrt-Schmidt, the second and third ones are bounded (the support 
of x' is near the origin). Hence the trace norm is polynomially bounded in A, and the exponential 
decay of p2 will do the rest. 

The proof of (|3.22p does not contain any new ingredients. The ideas are the same: introduce 
(13.121) in (13. 4p . integrate by parts with respect to z, propagate exponential factors over resolvents, 
and use the fact that the distance between the support of x' and the support of $l, 1 — <&l or $' L 
is of order L. □ 

Corollary 3.5. The linear response contribution to the current admits the thermodynamic limit 
and 

lim = -w? / X 'MTr w ( e ^ H y(Q)e- isH [p(H),n ± ])d S . (3.25) 

Proof. This is a direct consequence of Corollary 13. 21 Proposition l3.31 and the Lebesgue dominated 
convergence theorem. □ 

4 The Adiabatic Limit 

From now on the leads are semiinfinite; the thermodynamic limit has been taken. Hence in (|2.6p 
one has to interpret L as infinite. The adiabatic limit rj — > in the formula (|3.25p will in fact be 
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an Abel limit once we can show that the limit lim e lsII V(Q)e lsIi exists, at least for the strong 

s — >oo 

topology and on a certain subspace of M[H]. In view of the definition of V(Q) (see (|2.7p ) we are 
reduced to studying the limits lim e* sH n±e _IsH . 

s — >oc 

o 

The idea is to consider the Hamiltonian H as the perturbation of a decoupled Hamiltonian H a 

o 

that commutes with the projections H±. We define a decoupled Hilbert space Ti a := (II_7i) © 
(lIo7i) © (n + 7i). The decoupled Hamiltonian is simply H with two extra Dirichlet boundary 

o 

conditions defined by the previous splitting. Let us denote it by H a . We assume that the internal 
boundary defined by the condition x\ — ±a is smooth enough; it is so in the cylinder case. It is 
important to note that: 

[H o ,n±]=0. (4.1) 
With these notations and definitions we now have: 

e isH n±e- isH = e isH e- is ^U ± e is ^e- isH . (4.2) 

Under our assumptions, the Hamiltonian H in (|4.3p is nonnegative and has no singular continuous 
spectrum. Let us further assume that H has no embedded eigenvalues, and only a finite number 
of discrete eigenvalues of finite multiplicity located below the essential spectrum (which can arise 
from the geometry we chose for our system X^, [16]). Then let us denote by E a the finite 
dimensional orthogonal projection corresponding to a discrete eigenvalue A a of H, by the 
projection corresponding to its absolutely continuous spectrum, and by Eq :— © E a = 1 — E^ 

a<N 

the finite dimensional projection on its discrete spectrum. We can write: 

Tr„ (e lsU V(Q)e-* sU [p(H),n±]) = Tr n (E oc e tsU V(Q) e - lsIi E oc [p(H),II±]) 
+ J2 Tl '« (E a e isIi V(Q)e- lsII E 00 [p(H),n±]) 

+ J2 Tr « {E oc e^V{Q)e- lsli E a [p(H),n±]) 

a<N 

+ ^■n{E a e lsU V(Q)e' lsIi E fj [p(U),U ± }). (4.3) 

a<N,P<N 

Let us show that the last term does not contribute to the current after the adiabatic limit. We 
have the identity: 

Tr n {E a e lsU V(Q)e- lsH E fj [p(H), n±]) (4.4) 
= e^- x '^Tm (V(Q)Ef3[p(H.),IL±)E a ) . 

Now for a = (3 we have 

E a [p(U),U ± ]E a = p(\ a )(E a n±E a -E a Il ± E a ) = 0. (4.5) 
For a ^ P we use the bound 

\Tr H (V(Q)Ep[p(n),Il±]E a )\ < \\V(Q)\\ ||[p(H),II ± ]|| Bl 
and the fact that the s-integral 

V [ *'M e isu ds = (T X ')(u/v) ^ n 
J-oc 

for any ui ^ 0, as the Fourier transform of a L 1 -function on R. 
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Now let us study the second term in (|4.2p and prove that it will also disappear after the 
adiabatic limit. We have the identity: 

Tr n (E a e lsU V(Q)e- isU E 00 [p(H), n±]) (4.6) 
= Tr n (v(Q)e-^ H - x ^E OQ [p(H), U±] E« 



(4.7) 



and observe that (by interchanging the trace with the s-integral) we have 

Tr w jnQ) (v f x'(vs)e- ls{H - Xa) d S ^ E^ [p(H),U ± ] E Q 
= Trn{v(Q) (Vx'XT 1 (H - A«, )£«,)) [p(H), n±] E a } . 
Due to our hypothesis on H and the definition of Eoa, we get: 

(H-AaJ^oo > (inf ((T ac (H)) - A a )£oo, 
and since T\' converges to zero at infinity we obtain: 

\\W){r l &-\*)E°°)\\w\= su p (fx'Wv) — > o. 

ft>inf((Tae(H))-A a J '^ u 

The third term in ()4.3|) can be treated in a similar way as the second one. Therefore only the 

o 

first term can give a contribution, and let us identify it. Denote by P QC (H) = II_ © n + the 

o 

projector on the absolutely continuous subspace of H. We note that the incoming wave operators 

o 

at — oo associated to the pair of Hamiltonians (H, H) exist and are complete. This can be shown 
in a number of different of ways, but here we choose to invoke the invariance principle and the 
Kato-Roscnblum theorem. Indeed, the function — p is admissible (see Thm. XI. 23 [RS III]), and 
we have: 

n+(H, H) = n+(-p(H), -p(H)) = n+(-p(H), -p(H) - Ap), (4.8) 

o 

where the operator Ap := p(H) — p(H) is trace class (more details will be given in the next 
section). Thus: 



n+:=s- lim e »«M e -"Hp ac (H), Ran(f2+) = (4.9) 

S— * — OO 

Using JO), (Bill) and the fact that V(Q) = V{Q)P ac (H), we obtain: 
lim Tx n (E 00 e lsIi V(Q)e- tsIl E 00 [p(H),n±]) 

s — *■— OO 

= lim Ti n (E 00 e isI1 e- is ' k V(Q)e ls ' k e- lsU E 00 [p(H),n ± } 

= Tr„ (E^fl+ViQ^XE^ [p(H), n±]) . (4.10) 

Using (|4.10p and (|4. 3[) in (|3.25|) we can write a formula for the adiabatic limit of the linear response 
current at t = as: 

J± ■= lim j v = iTr n (E^n+ViQ^E^ [p(H),U±]) (4.11) 

•q— »0 

= v-iTr E „H (o+n_n; [p(H),n±]) + v+m Eoa H (o+n+n; [p(H),n±]) . 

Recalling that the decoupled Hamiltonian commutes with the projections II ± (see (J2H])), we can 
write: 



[p(H),n±] = 



(p(H)-p(H)),n± 
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Then we have: 

i± = v-iTr EeBn (n+n_n* + [A P ,n±]) +v + iTi ExH (n + u + n* + [a p ,il±]) . (4.12) 

At this point we can show that if v+ = V—, then the current is zero. Indeed, since f2+(II+ + 
n_)J7^ = 1 on EooTL, and because Tr-^([Ap, 11+]) = we can write: 

j+ = {v- - v+)iTr Eoon (0+n_n; [Ap,n+]) - iTr EppH ([A P ,IL + }). (4.13) 

But Tr EppH ([A/9,n+]) = Tr(£ pp (H) [p(H),II+]) = 0, therefore we obtain: 

3+ = («- - v + )iTr Eoon (ft + n_ft+ [Ap,II+]) , (4.14) 

which is a close analog of the formula given by NESS type approaches (see [I] and references 
therein). A remark is in order here. In NESS type approaches the reservoirs (leads) are suddenly 
coupled at t = and the current is computed in the limit t — > oo. In general, the current has 
an oscillatory component given by a term similar to the last term in (|4.3[) |19j . so it does not 
have a definite limit at t — > oo. What it is actually computed by an averaging procedure, is the 
steady component of the current. As we have seen above, in the adiabatic switching formalism the 
current itself has a limit as rj — > 0. This fact confirms (in our particular setting) the basic feature 
of the adiabatic switching, namely that it "kills" oscillatory terms of various physical quantities. 



5 The Landauer-Biittiker formula 

o 

We will now compute the trace appearing in equation (|4.14|) using a spectral representation of H. 
By doing this we will obtain the L-B formula. 

Let us enumerate a few technical results needed in order to justify the computations which will 
follow next. In the next subsection we will only give the proof of (|5.ip . since the other points are 
just well known facts which can be found in any standard book treating the stationary scattering 
theory. 

Proposition 5.1. Consider some n > 1/2. 

o 

1. The operator Ap = p(H) — p(H) is trace class, and if (x) := y/l + x\ we have: 

0"(Ap)O n eB(tt). (5.1) 



2. Denote by Eq > the infimum of the essential spectrum of H. Then er ac (H) = \Eq,oo). 
Moreover, the spectrum o/H in [Eq, oo) is absolutely continuous, with finitely many (possibly 
none) eigenvalues. The set of thresholds T (points where we do not have a Mourre estimate, 
see \21$ ) is discrete and contains infinitely many points. Ifn is large enough then the mapping 

A : [E a , oo)\Th B(W), (5.2) 
A(X) := lim (.)- ,l (H-A- J£ )- 1 (.)-" 

is continuously differentiable. 

3. For n large enough, the mapping 

A p :[0,p(Eq)}\p(T)^M(H), (5.3) 



A p (t) := lim (.)- n (p(H)-t- l e) 

e\0 



-1 , 



is continuously differentiable. 
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4- We can choose a family of generalized eigen] 'unctions for the absolutely continuous part of the 

o 

operator H, which in the case of our cylinder-like semiinfinite leads can be explicitly written 
down. The multiplicity of the absolutely continuous spectrum is finite when the energy is 
restricted to compacts, but it is not constant and increases with the energy. A given channel 
will be indexed by {a, a), where a e {±1} shows the left/right lead, and a indexes all other 
quantum numbers at a given energy. Hence for n > 1/2 we have : 

o(a,±) o (os,±) oK±) 

(■)-% eU, Hv E '"="E%, E ', Ee[E ,^)\T. (5.4) 

5. With this choice, we can define the corresponding generalized eigenfunctions of H either with 
the help of the limiting absorption principle: 

{■)-> { e ±] = (r^E^ mp(e)) {(r^p{r}{(r n ^E ±} ) , (5.5) 

Ee [Eo,oo)\T, 

or as solutions of the Lippmann-Schwinger equation (with the usual abuse of notation): 

<p%<V = - ( P (H) - p(E) l + y\A P )^ ± \ (5.6) 

Ee [E ,oo)\T. 

If n is large enough, then the map 

[E ,^)\T3E^(-)- n ^ ±] £H (5-7) 

is continuously differentiable. 

Let us now use this proposition in order to compute the current. The main idea is to compute 
the integral kernels of the two trace class operators in (|4.12[) , and then to compute the trace as the 
integral of the diagonal values of their kernels. We use for this the generalized eigenvalues of H. 
Let us choose two energies E, E' e er ac (H) \7~. Denote by (a, ±) and (a', ±) the quantum numbers 
describing the spectral multiplicity of H at E and E' . Then in the spectral representation of H, 
an operator like £' oo r2 + II_r2^(A/9)II + will have the integral kernel (here er, a' e {+, — } indicate 
the leads): 

l(E,a,o-,E',a',o-') = 5 CT , + ((A /J )n + ^'' <T ' ) , ^' +) }- (5.8) 
In deriving this equation we formally used the " intertwining" property 

n + tp E = <p E (5.9) 

o 

since the wave operators are unitary between the absolutely continuous subspaces of H and H. 
The scalar product is in fact a duality bracket between weighted spaces, and we used (|5.2p . This 
integral kernel is jointly continuous in its energy variables outside the set of thresholds, due to 
(|5.7|) . Thus when we compute the trace of i? oo r2 + n_0^(A ( o)n + , we may write: 

Tr-E 00 f2 + n_f2+(Ap)II + = lim y^l(E,a,o;E,a,o)dE, (5.10) 

S-.» ac (H)\Tj s a a 

where S denotes compact sets included in <r ac (H) \ T, and the limit means that the Lebesgue 
measure of er ac (H) \ (T U S) goes to zero. 

Now after undoing the commutator in (|4.14|) . we compute the traces in the way we described 
above, keeping in mind that the energy integral has to be understood as a limit avoiding the 
thresholds. Thus we obtain: 



16 



-« + )([Ap,n + ]^'- ) ,^'- ) )} 

a 

/>OC 

= 2(v--v + ) Y.{^ A P U + ( PE~\'P { E~ ) )}dE. (5.11) 

Using the Lippmann-Schwinger equation (|5.6p we obtain 

ApU+cp^ = -ApIL + (p(h) - p{E) - i + )- 1 (A P )^-\ 

o 

and we insert this in (|5.11[) . Using the generalized eigenfunctions of H, the limiting absorption 
principle, the various regularity properties listed in Proposition 15. 1[ and the Sokhotskii-Plcmelj 
formula, we obtain: 



3 ( n+(p(H) - p(E) - l Q + )- 1 (Ap)^-\(Ap)^ 



,(«■-) 



(5.12) 



o («',+) 



Thus the current reads as: 

j+ = 2n{v+ 



{Ap)^ E ',<p E 



dE. 



(5.13) 



We now have to relate the above integrand with the scattering matrix. We know that the S matrix 

o 

commutes with H, and so does the T matrix defined as ^(1 — S). In the spectral representation 

of H induced by {ip E } a ,a, the T operator is a direct integral with a fiber which is a finite 
dimensional matrix. Now using the correspondence principle (see (|4.8jl ). and formula (4) ± , page 
233 in [H] we can write (see also Thm. XI.42 in (2SJ): 



1 a.cr-a' ,a 



{e) = 7(e)\ { Ap) ^ E ' 



The Landauer-Biittiker formula becomes: 



j+ = 2tt(v + -v.) Yl \T a ',+-, a ,-(E)\ 2 dE. 

•I e q , 



(5.14) 



(5.15) 



We can now state the main result of our paper: 

Theorem 5.2. Assume that the system has N leads. Then the charge variation of lead f = 1 
defined after the thermodynamic and adiabatic limits, in the linear response regime, reads as: 



N 
/=2 



(5.16) 



In terms of transmittance between leads we have: 

N 

Jl 



POO 

2irJ2( v i-Vf) / ^'(£)|T a ,, w (£)| 2 <iE. 

f=2 JE ° a,a' 



(5.17) 
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5.1 A continuity equation for the current 

If we go back to (|2.35[) and instead of H± we use some smaller projections II± which have the 
property that n±II± = IT ± (but still corresponding to an infinitely long portion of a lead), then 
we can again define a current and obtain an expression as in (|2.44[) , where Tl± are replaced by II± . 
The thermodynamic and adiabatic limits can be performed in the same manner, and we would 
obtain a formula very similar to (15.161) : 

N 

ji = E( v / -vi^E^n (n+n/n^ro.nx]) . (5.18) 

o 

Note that H is still defined with respect to the old decomposition imposed by IPs, and NOT by 

ill. That is H does NOT commute with fix . Another important observation is that [p(H),fli] is 
already trace class even before inserting something commuting with IT, as we can see from the 
proof of Proposition 13. 31 If wc subtract ji from ji we obtain (denote by ~k\ := Hi — Hi): 

N 

h-h = Y,( v f- Vl ) iTtB -« (« + IW[ p (H),^]) . (5.19) 

/=2 

Now 7Ti corresponds to the multiplication with a function with compact support. Using again the 
methods of Proposition ^. 3[ one can show that p(H.)wi is trace-class (one can write it as a product 
of two Hilbert-Schmidt operators). Then 

TWiCfi+n/fi+bCH),*-!]) =Tr(0 + n / fi; /9 (H) 7 r 1 )-Tr(0 + n / 0;vri / 9(H)). (5.20) 

Let us show the operator /3(H)Q+ri/fi5_7ri is also trace class: note first that [O+II/Ql, /o(H)] = 

o 

(this is due to the intertwining property of wave operators and due to [H, II/] = 0); second, we can 
put near w± enough resolvents in order to obtain the product of two Hilbert-Schmidt operators. 
Then let us prove the trace cyclicity: 

Tr (ft+iP/f^TTipQH)) = Tr (p(H)Q + Il f n* + TTi) . (5.21) 

Indeed, take ~Ki with compact support such that tt2^i = 7Ti. The usual trace cyclicity gives: 

Tr (n + U f n* + nip(U)) = Tr (ft+IL/fi^TriTra^H)) = Tr (^(HjO+IL/fi^Tri) . (5.22) 

Now take back tt2 to the right and (|5.21[) is proved. Thus we can write: 

Tr (Q + U f n* + TTip(U)) = Tr (/>(H)0 + II/f2* vri) = Tr (n + U f n* + p(H)ni) (5.23) 

where in the second equality we used that [0 + ri/0^_, /o(H)] = 0. Use this in (|5.20[) and obtain 

h = 3i, 

which shows that the current is the same no matter where we measure it on the lead. 

5.2 Sketching the proof of (ED) 

Proof. The main idea is to use the geometric perturbation theory, in the same spirit as in Propo- 
sition 13.31 With the functions introduced in (|3.6|) and (|3.7J) , choose L large enough such that 

H$l = H$i. This means that the "sample" is contained in the domain where = 1. Then 
keep L fixed. Now for z £ C \ K we define the parametrix: 

S(z) := $ L (H - z)- 1 ^ + $ L (H - zy^L. (5.24) 
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We have that the range of S(z) is mapped into the domain of H and: 



(H-z)S(z) =: 1 + T(z), 

T(z) = [-2(V*i)V - (A$i)](H - zY^l 



(5.25) 



+ [-2(V*i)V - (A$ L )](H - z)- 1 ^. 



Similarly as in (|3.12p we can write: 



(H - z) 



-i 



S{z)-(H-z)- 1 T(z). 



(5.26) 



Thus: 



(H-z) 



(H - z)- 1 = $ L (H - z)- 1 $ L + ($£ - 1)(H - z)- 1 ^ 



+ (H-z)- 1 ($ L -l)-(H-z)- 1 T(z). 



(5.27) 



Then using the Cauchy integral representation as in (|3.15p , we see that the operator Ap can be 
expressed as an integral where the integrand consists of several terms, and we can write each of 
them as a product of two Hilbert-Schmidt operators, due to the support properties of our cut-off 
functions and the exponential decay of various integral kernels. In order to prove (|5.ip we need 
to propagate some exponential decay near both polynomials (■)", and this can be done as in 
Proposition 13. 31 We do not give further details. □ 
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